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Recently we reported that the BH3-only proteins Bim and Noxa bind tightly but transiently to the 
BH3-binding groove of Bak to initiate Bak homo-oligomerization. However, it is unclear how such tight 
binding can induce Bak homo-oligomerization. Here we report the ligand-induced Bak conformational 
changes observed in 3D models of Noxa - Bak and Bim -Bak refined by molecular dynamics simulations. In 
particular, upon binding to the BH3-binding groove, Bim and Noxa induce a large conformational change of 
the loop between helices 1 and 2 and in turn partially expose a remote groove between helices 1 and 6 in Bak. 
These observations, coupled with the reported experimental data, suggest formation of a pore-forming Bak 
octamer, in which the BH3-binding groove is at the interface on one side of each monomer and the groove 
between helices 1 and 6 is at the interface on the opposite side, initiated by ligand binding to the 
BH3-binding groove. 



Paralogs of the B-Cell lymphoma 2 (BCL-2) gene express proteins that regulate the cellular life-or-death 
switch by controlling mitochondrial outer membrane permeabilization (MOMP) 1 " 5 . Sharing up to four 
conserved BCL-2-homology (BH) domains (viz., BH1, BH2, BH3, and BH4) 1,6 , these proteins can be 
divided into anti- and pro-apoptotic subgroups. The anti-apoptotic members such as Bcl-2 contain four BH 
domains and inhibit MOMP. The pro-apoptotic members, all of which promote MOMP, can be further divided 
based on structural and functional differences into effectors and BH3-only proteins. Effectors, including Bax (Bcl- 
2-associated x protein) and Bak (Bcl-2 antagonist killer), contain BH1, BH2 and BH3 domains and are capable 
of directly causing MOMP. In contrast, the BH3-only proteins exhibit homology to Bcl-2 only in their BH3 
domains, produce MOMP indirectly, and can be subcategorized into activators and sensitizers 7 . The activators 
such as Bim (Bcl-2-interacting mediator of cell death) appear to bind Bax or Bak to trigger homo-oligomeriza- 
tion 8,9 . The sensitizers such as Bad (Bcl-2 antagonist of cell death) are thought to indirectly promote MOMP by 
binding and neutralizing the anti-apoptotic members 1 . 

Protein-protein interactions involving various Bcl-2 paralogs play a critical role in their function. The binding 
of BH3-only proteins to the anti-apoptotic proteins involves a surface groove called the BH3-binding groove as 
revealed by crystal or nuclear magnetic resonance (NMR) structures of anti-apoptotic protein complexes. The 
effectors also have a groove akin to the BH3-binding groove of the anti-apoptotic proteins. Whether this groove in 
Bak or Bax binds BH3 domains of direct activators and contributes to MOMP has previously been unclear. 

Although the BH3-only protein Noxa (Latin for damage) was widely considered to be a sensitizer 7,1012 , it has 
been recently reported that Noxa can function as a direct activator of Bak and Bax 13,14 . We have further reported 
that, unlike the sensitizer Bad, Noxa and Bim bind tightly (K D values of 20-30 nM) but transiently to the BH3- 
binding groove of Bak, thereby triggering the formation of a Bak homo-oligomer of —200 kDa, MOMP, and 
killing of intact cells 14 . 

Despite the fact that Bak and Bax play pivotal roles in apoptotic responses and new information about their 
regulation by BH3-only proteins has become available 1315 , there has been no structural information on interac- 
tions of BH3-only proteins with the BH3-binding groove of Bak or Bax except for a brief report of a computational 
three-dimensional (3D) model of the Noxa - Bak complex 14 . It has been unclear how a high-affinity activator can 
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bind Bak transiently and trigger Bak ftomo-oligomerization, as one 
would expect hetero-oligomerization to be induced by binding of 
Noxa to Bak with a K D value of 24 nM. 

In this report, we describe the ligand- induced Bak conformational 
changes observed in computational 3D models of Noxa "Bak and 
BinrBak that were refined by near-microsecond-scale molecular 
dynamics simulations. In the context of these complex models and 
the Bak conformational changes, we discuss the nature of interac- 
tions of BH3-only proteins with the Bak BH3-binding groove and a 
possible mechanism for Bak homo-oligomerization induced by tight 
but transient binding of activator BH3-only proteins. 

Results 

The Noxa - Bak Complex Model. Unless specified otherwise, human 
Bcl-2 paralogs were used in our studies and the Bak helical 
nomenclature of reference 16 is used in this report. Our homology 
modeling using the SWISS-MODEL program 17 predicted that 
residues 19-45 of Noxa 18 adopt a helical conformation, a finding 
that has subsequently been confirmed in an unpublished crystal 
structure of the helical Noxa binding at the BH3 -binding groove of 
the anti-apoptotic paralog Bcl-2-related protein Al (BFL-1) 19,20 
released after this study (Protein Data Bank ID: 3MQP). Based on 
our initial site-directed mutagenesis data involving residues at the 
Bak BH3-binding groove 14 and the structural similarity between the 
apo Bak crystal structure and the NMR structure of mouse anti- 
apoptotic protein Mcl-1 (myeloid cell leukemia- 1) whose BH3- 
binding groove is occupied by a mouse Noxa BH3 domain 
(mNoxa) 21 , we conjectured that the helical Noxa would dock at the 
BH3-binding groove of Bak, even though the groove is occluded by 
R88 Bak and Y89 Bak in all the reported apo Bak crystal structures 
(Figure la) 22,23 . We then generated the Noxa - Bak complex by 
manually docking the helical Noxa atop the BH3 -binding groove 
of a Bak conformer (residues 21-183) that was taken from a crystal 
structure of an apo Bak homo-dimer with the C-terminal trans- 
membrane domain excised 22 . The manual docking was performed 
in such a way that L29 Noxa was placed as close to I114 Bak and L118 Bak 
as possible and that L36 Noxa was as close to V129 Bak and I85 Bak as 
possible. This model was refined first by energy minimization, 
which moved Noxa closer to the BH3-binding groove and pushed 
R88 Bak out of the groove. The resulting complex was further refined 
by multiple molecular dynamics simulations using a published 
protocol 24 . Seventy-two independent 10-ns simulations were 
carried out, thus offering collectively 0.72-us sampling for the 
Noxa - Bak conformations. Although Y89 Bak was still within the 
BH3-binding groove and underneath Noxa in the energy- 
minimized Noxa "Bak model, as described below, the model refined 
by the simulations has Y89 Bak pushed out of the groove and the BH3- 
binding groove expanded to accommodate Noxa fully (Figures lb 
and lc). To provide the information on the protein dynamics and the 
solvent effect on complex stability as described below, we release 
coordinates of average and energy-minimized average confor- 
mations of Noxa - Bak (datasets SI and S2, respectively) as well as 
coordinates of a water-containing instantaneous conformation 
of Noxa - Bak that has the smallest root mean square deviation 
(RMSD) to the average conformation (dataset S3). 

Computational Support for the Noxa - Bak Model. As a control 
study, we performed 11 10-ns simulations of the first model of 
the NMR structure of mouse Mcl-1 in complex with mNoxa 
(mNoxa - mMcl-1) 21 using the same simulation protocol as the one 
for Noxa'Bak. We also carried out 11 10-ns simulations of the first 
model of the NMR structure of mouse Mcl-1 in complex with mouse 
Puma (p53 up-regulated modulator of apoptosis; mPuma - mMcl- 
l) 21 . These simulations rendered two average conformers with 
alpha carbon root mean square deviations (ACRMSDs) of 1.5 A 
for mNoxa - mMcl-l and 1.8 A for mPuma - mMcl-l relative to the 



respective average NMR structures (Table 1). Each of the two 
computational average conformers was obtained by averaging all 
11,000 trajectories saved at 1.0-ps intervals during the last 1-ns 
period of the 11 10-ns simulations. Reassuringly, these ACRMSDs 
are less than the corresponding ones of the 20 NMR models of the 
mNoxa - mMcl-l and the 20NMRmodels of mPuma - mMcl-l (2.5 A 
for mNoxa - mMcl-l and 2.0 A for mPuma - mMcl-l). These results 
suggest that the 0.72-us molecular dynamics simulation protocol is 
applicable to Noxa - Bak, because the sizes of mNoxa - mMcl-l (189 
residues) and mPuma - mMcl-l (189 residues) are comparable to that 
of Noxa - Bak (190 residues). 

Because the Bak BH3-binding groove is mainly hydrophobic, heli- 
ces such as polyalanine can bind to the groove with favorable inter- 
molecular interaction energies. To be biologically relevant, the 
intermolecular interaction energy of Noxa with the Bak groove 
should intuitively be lower than that of polyalanine. To confirm this, 
we calculated intermolecular interaction energies of the Noxa - Bak 
complex and a reference complex of A 126 K 27- Bak that was generated 
by mutating residues 19-45 of Noxa in the Noxa - Bak complex model 
to A 19 - 44 K 45 . The mutation to A 19 - 44 K 45 rather than A 19 45 was to keep 
the net charge of the mutant peptide identical to that of Noxa. To 
obtain lower intermolecular interaction energy for the mutant com- 
plex, we refined the complex by running 1 1 10-ns simulations using 
the same simulation protocol for Noxa - Bak. Reassuringly again, 
the intermolecular interaction energy and the interaction energy 
per residue of the reference complex are —110.0 kcal/mol and 
— 4.1 kcal/mol, respectively (Table 1), while the corresponding 
energies for Noxa - Bak are —132.0 and —4.9 kcal/mol, respectively. 
These energies suggest that the Noxa • Bak model is worthy of further 
experimental testing. 

Experimental Support for the Noxa - Bak Model. Our sequence 
analysis identified L29 and D34 as two highly conserved residues 
in the Noxa BH3 domain and L29, F32, and L36 as the hydro- 
phobic residues in this domain 14 . The double (L29A and D34A) 
and triple (L29E, F32E, and L36E) mutants of Noxa showed 
diminished binding to Bak 14 , suggesting that the identified residues 
are crucial for Bak binding. Our sequence analysis also identified 
R127 as a highly conserved residue in the BH3-binding grooves of 
Bak, Bax, Bcl-2, Bcl-x L (Bcl-2 X protein), and Mcl-1. The Bak R127A 
mutant exhibited markedly reduced binding to Noxa 14 , indicating 
that BH3-binding groove is involved in the binding of Bak with 
Noxa. Acquired prior to the Noxa - Bak complex prediction, these 
mutation data are consistent with the Noxa - Bak model in view of 
the pair-wise intermolecular interactions between Noxa and Bak 
as detailed in Table 2 and Figures Id- If. In addition, the model 
suggested that F93, 1114, and G126 of Bak directly interact 
with Noxa (Table 2 and Figures ld-lf), and subsequent single 
mutations of these residues (F93E; I114E; G126S) markedly 
reduced the binding of Bak with Noxa 14 . Although NMR or X-ray 
crystal structure of Noxa - Bak is required to further validate the 
model, the mutagenesis data suggest that the computational model 
can be used with care for obtaining insight into Noxa-induced 
Bak conformational changes that are crucial to the Bak homo- 
oligomerization. 

Bak Conformational Changes Induced by Noxa. Expanding the 
Bak BH3-Binding Groove. Compared to the two reported crystal 
structures of apo Bak 22,23 , one noticeable change of the Bak in 
complex with Noxa is main-chain and side-chain conformational 
changes of R88 Bak and Y89 Bak (Figure 1). In the apo Bak crystal 
structures, Y89 is in a network of 10 aromatic residues that 
interact with one another via pi-pi interactions (Figure 2a). This 
aromatic network, which includes Y89, Y108, Y110, F93, Fill, 
F119, F134, F150, F157, and F161, presumably provides structural 
stability to Bak in water. Hence one would expect that, upon Noxa 
binding, R88 Bak might move away from the BH3 -binding groove and 
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Figure 1 | Bak conformational changes induced by binding of Noxa to the BH3-binding groove, (a) The BH3-binding groove of the unbound Bak 
showing the blockage of the groove by R88 bak and Y89 bak . (b) The BH3-binding groove of the bound Bak showing the vacated groove by R88 bak and Y89 b * k 
revealing two hydrophobic holes, (c) The helical Noxa anchored at the Bak BH3-binding groove showing complementarity between Noxa and the groove, 
(d) The Bak residues that interact with Noxa. (e) The Noxa residues that interact with Bak or solvent, (f) Intermolecular interactions between Noxa and 
the Bak BH3-binding groove and between F38 No "* and solvent molecules (not shown) . (g) The unbound Bak with a water molecule hydrogen-bonding to 
Bak residues underneath the turn between two helices, (h) The bound Bak with a water-filled cavity underneath the turn between two helices, (i) The 
bound Bak with a 7-water hydrogen-bond network involving residues of Bak and Noxa in the cavity underneath the turn between two helices. The 
unbound Bak is taken from the crystal structure of Protein Data Bank code of 2IMT, and the bound Bak is taken from the simulation-refined Noxa - Bak 
model. For clarity hydrogen atoms are not displayed except for the nonpolar hydrogen atoms of Gly. Atoms of oxygen, nitrogen, sulfur, hydrogen, carbon 
of Bak, and carbon of Noxa are colored in red, blue, yellow, magenta, gray, and green, respectively. Dashed lines denote hydrogen bonds. 



Y89 Bak might remain in the groove to maintain the integrity of the 
aromatic core. Surprisingly, Noxa pushed Y89 Bak out of the groove 
through main-chain and side-chain conformational changes of Bak. 
These changes are accompanied by expansion of the groove and 
formation of two hydrophobic holes in the groove (Figure lb), one 
comprised of I81 Bak , I85 Bak , and V129 Bak , and the other made of 
F93 Bak , I114 Bak , L118 Bak , and A130 Bak . Consequently, Noxa snugs 
down to the groove with its L29 Noxa and L36 Noxa plugged into the 
two holes. This groove expansion also creates a cavity (Figure lh) 
allowing seven water molecules to form a hydrogen-bond network 



underneath the turn that connects helices a2 and oc3 of Bak 
(Figure li; dataset S3), whereas only one water molecule is present 
underneath the turn (Figure lg) as revealed by the high-resolution 
crystal structure of the apo Bak 22 . The 7-water hydrogen-bond 
network, which involves residues from both Noxa and Bak 
(Q28 Noxa , G82 Bak , I85 Bak , N86 Bak , Y89 Bak , and G133 Bak ; Figure li), 
contributes to the affinity of Noxa for Bak. 

Expansion of the Aromatic Network ofNoxa'Bak. As shown by the 
energy-minimized average conformation of the Noxa - Bak model 
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Table 1 | Interaction Energies and Root Mean Square Deviations of Bak or mMcl-1 Complexes 



Complex 

Noxa-Bak 

A i-26 K 27. Bak 

mNoxa-mMcl-1 
mPuma-mMcl-1 



E (kcal/mol) 

-132 
-110° 
-176° 
-152= 



E / BH3 residue 
(kcal/mol) 

-4.9 
-4.1 ° 
-6.5° 
-5.8° 



CRMSD (A) 

1.2 b 
1.8 C 
1.5 d 
1.8 d 



ACRMSD"™ 5 
(No. of conformers) (A) 

2.7(5,100)" 
2.4 (2,000) F 
2.8 (1 1,000)° 
2.6 (1 1,000)° 



ACRMSD NMR 
(No. of conformers) (A) 

NA 

NA 
2.5 (20) h 
2.0 (20) h 



Q E: the intermolecular interaction energy and the intermolecular interaction energy per residue of BH3-only protein were obtained using the same method as the ones for Noxa-Bak except that 1 1 10-ns-long 
simulations were performed. The initial structure of mNoxa-mMcl-1 ormPuma-mMcl-1 used for the simulations was taken from the first NMR model of Protein Data Bank code of 2ROD or 2ROC, respectively. 
The initial structure of A'~ 26 K 27 -Bak was obtained from mutation from the simulation-refined Noxa-Bak model. 

b CRMSD: alpha carbon root mean square deviation between the average structure of 5,1 00 similar conformers identified by cluster analysis from the second-round simulation of Noxa-Bak and the energy- 
minimized average structure of 271 similar conformers identified by cluster analysis from the first-round simulation of Noxa-Bak, where the average structure of the 5, 1 00 conformers was not subjected to 
any energy minimization. 

C CRMSD: alpha carbon root mean square deviation between the average structure of 2,000 similar conformers identified by cluster analysis from the simulations of A'~ 26 K 27 - Bak and the energy-minimized 
structure of A'~ 26 K 27 -Bak derived from mutations of the energy-minimized average structure of the 5,100 Noxa-Bak conformers. 

d CRMSD: alpha carbon root mean square deviation between the average structure of 1 1 ,000 similar conformers identified by cluster analysis from the simulation of mNoxaA-Mcl-1 (or Puma-Mcl-1 ) and the 
average structure of the 20 NMR models of mNoxaA-Mcl-1 (or Puma-Mcl-1 ), where the two average structures were not subjected to any energy minimization and the cluster analysis identified only one 
cluster of conformers. 

e ACRMSD MMDS : average alpha carbon root mean square deviation among the 5,100 similar conformers of Noxa-Bak. 
'ACRMSD^ 05 : average alpha carbon root mean square deviation among the 2,000 similar conformers of A'~ 26 K 27 -Bak. 
9 ACRhASD MMDS : average alpha carbon root mean square deviation among the 1 1 ,000 similar conformers of mNoxa-Mcl-1 (or Puma-Mcl-1). 
h ACRMSD NMR : average alpha carbon root mean square deviation among the 20 NMR conformers of mNoxa-Mcl-1 (or Puma-Mcl-1 ). 



(Figure 2b and dataset S2), F32 Noxa occupies the groove region where 
Y89 Bak resides in the apo Bak conformation. Consequently, the 10- 
residue aromatic network is expanded to an 1 1 -residue network after 
the groove is occupied by Noxa. Furthermore, as revealed by the 
average conformation of the Noxa - Bak model (dataset SI), 8 of 
the 11 aromatic residues (F32 Noxa , Y89 Ba \ Y108 Ba \ F93 Bak , Flll Bak , 
F134 Bak , F150 Bak , and F157 Bak ) showed ring contraction due to ring 
spin with or without side-chain motion. In particular, the side chains 
of F32 Noxa , Y89 Bak , F93 Bak , and F134 Bak are markedly contracted rela- 
tive to the side chain of Y110 Bak (Figure 2c). These observations 
suggest that the high affinity of the Noxa - Bak complex stems from 
both enthalpy and configurational entropy 25 . 

Exposing the Groove between Helices al and a6. The most noticeable 
Bak conformational change induced by Noxa occurs in the loop 
between helices al and oc2 (residues 51-61; Figures 3a and 3b), which 
in turn partially exposes the groove between helices al and a6 that is 
located on the opposite side of Bak from the BH3 -binding groove 
(Figure 3b). This is caused by the Bak BH3 -binding groove expan- 
sion, which pushes helix a5 slightly away from helix a3 (Figure 3c). 
In the apo Bak crystal structure, because the BH3-binding groove is 
contracted, there is room between helices al and a5 to accommodate 



residues 51-61 (Figure 3c). In the bound Bak computational model 
(Figure 3b), because the BH3-binding groove is expanded, there is no 
room between helices al and a5 for residues 51-61, which leads to 
the large conformational change of residues 51-61 and in turn par- 
tially exposes the groove between al and a6 (Figure 3c). As discussed 
below, partial exposure of this groove might be one of the key steps in 
Bak homo-oligomerization. 

The BinrBak Model and Bim-Induced Bak Conformational 
Changes. The computational 3D model of the BinrBak complex 
was subsequently generated by mutating Noxa (residues 19-45) of 
the Noxa - Bak model to Bim (residues 143-166) followed by 
refinement with energy minimization and then 100 10-ns mole- 
cular dynamics simulations using the same simulation protocol as 
the Noxa - Bak model. Coordinates of average and energy-minimized 
average conformations of the Bim - Bak model are provided in 
datasets S4 and S5, respectively; coordinates of a water-containing 
instantaneous conformation of Bim - Bak that has the smallest RMSD 
to the average conformation are provided in dataset S6. As apparent 
from datasets S4-S6, induced by binding of Bim to the Bak BH3- 
binding groove, the Bim - Bak model has almost the same Bak 
conformational changes as those of the Noxa - Bak model, except 
that the Bim-induced conformational change of the loop between 
helices al and a2 exposes more of the groove between al and a6 in 
two Boltzmann-weighted minor conformations than the Noxa- 
induced loop conformational change (datasets S7 and S8). 

Discussion 

Computational models of Noxa - Bak and Bim - Bak. To date there 
has been no experimental structure of Bak or Bax with its BH3- 
binding groove occupied by a BH3-only protein. An NMR struc- 
ture has been reported for Bax with its C-terminal transmembrane 
domain bound at the BH3-binding groove and a conformationally 
constrained peptide of the Bim BH3 domain docked at a groove 
between two helices that correspond to helices al and a6 of Bak 15 . 
Derived from manual docking and extensive refinement with near- 
microsecond-scale molecular dynamics simulations, the computa- 
tional model of Noxa bound in the BH3-binding groove of Bak is 
supported by the computational self-consistency tests in terms of 
conformational sampling and intermolecular interaction energy 
and by site-directed mutagenesis studies. This model is also sup- 
ported by an unpublished crystal structure of Noxa in complex with 
BFL-1 (Protein Data Bank ID: 3MQP; released after our model pre- 
diction); the alpha-carbon and heavy-atom RMSDs of Noxa (resi- 
dues 19-40) between the crystal structure and the energy-minimized 



Table 


2 1 Resid 


ues Involving Favorable 


Intermol 


ecular Interactions 


between Noxa 


and Bak 








IDNoxa" - IDBak 


Interaction 


IDNoxa" - IDr^ 


Interaction 


C25- 


-M96 


van der Waals 


F32- 


185 


van der Waals 


C25- 


-Yl 10 


van der Waals 


F32- 


Y89 


pi-pi 


C25- 


-11 14 


van der Waals 


G33- 


-LI 18 


van der Waals 


A26- 


-1114 


van der Waals 


G33- 


-R127 


van der Waals 


A26- 


-SI 17 


van der Waals 


G33- 


-A130 


van der Waals 


Q28- 


-E92 


van der Waals 


D34- 


-R127 


ionic 


Q28- 


-M96 


van der Waals 


K35- 


-Y89 


cation-pi 


L29- 


M96 


van der Waals 


L36- 


185 


van der Waals 


L29- 


Fl 1 1 


van der Waals 


L36- 


G126 


van der Waals 


L29- 


11 14 


van der Waals 


L36- 


V129 


van der Waals 


L29- 


Al 15 


van der Waals 


L36- 


A130 


van der Waals 


L29- 


LI 18 


van der Waals 


N37- 


-G126 


van der Waals 


L29- 


A130 


van der Waals 


N37- 


-G126 


hydrogen bond 


L29- 


L131 


van der Waals 


R39- 


-D84 


ionic 


L29- 


F134 


van der Waals 


Q40- 


-W125 


hydrogen bond 


R30- 


-LI 18 


van der Waals 


L43- 


181 


van der Waals 


F32- 


F93 


pi-pi 


L43- 


D84 


van der Waals 



SCIENTIFIC REPORTS | 2 : 257 | DOI: 1 0. 1 038/srep00257 



4 




Figure 2 | Expansion and side-chain dynamics of the aromatic core inside 
the Noxa-Bak complex, (a) The 10-residue aromatic core of the apo Bak 
crystal structure (Protein Data Bank code: 2IMT). (b) The 11-residue 
aromatic core of the Noxa - Bak model obtained with energy minimization 
after molecular dynamics simulations, (c) Close up view of the highly 
mobile aromatic residues of the Noxa - Bak model obtained without energy 
minimization after molecular dynamics simulations. 



average conformation of the model (dataset S2) are 0.9 and 2.3 A, 
respectively. Energetically, one may have a concern with the 
Noxa'Bak model in which F38 Noxil is exposed to solvent rather than 
being buried in the BH3-binding groove (Figures le and If). In the 
3MQP crystal structure at a resolution of 2.24 A, F38 Noxa is also 
exposed to solvent, despite a penalty of the solvation energy of the 
solvent-exposed phenyl group. While verification by X-ray crystal- 
lographic or NMR studies is necessary, the computational Noxa 1 Bak 
model and its derivative model of Bim ■ Bak offer insights as discussed 
below. 

Binding of BH3-Only Proteins to the BH3-Binding Groove. 
Binding of BH3-only proteins to the BH3-binding groove of anti- 
apoptotic proteins has been extensively studied. According to cur- 
rent literature the interactions with the anti-apoptotic BH3-binding 
groove include hydrophobic and ionic interactions 5,21,26-29 . In the 
Noxa - Bak model, Noxa interacts with hydrophobic residues of 
I81 Bak , I85 Bak , F93 Bak , V129 Bak , II 14 Bak , LI 18 Bak , and A130 Bak and ionic 
residues of D84 Bak , E92 Bak , and R127 B!lk at the BH3 -binding groove 
(Figures Id— If). These hydrophobic and ionic interactions are con- 
sistent with the interactions involving the anti-apoptotic BH3-bind- 
ing groove. 

Two other types of interactions also appear to contribute to the 
tight binding between Noxa and Bak. First, Noxa also interacts with 
Bak via pi-pi interactions involving the 1 1 -residue aromatic network 
of the complex. Second, cation-pi interactions also contribute. In the 
Noxa-Bak model, Y89 Bak , R88 Bak , K35 Noxa , and R39 Noxa are exposed to 
solvent. It is, therefore, difficult to discern from an instantaneous 
conformation or an energy-minimized average conformation whe- 
ther Y89 Bak forms cation-pi interactions with R88 Bak , K35 Noxa , and/or 
R39 Noxa . In the average conformation of the Noxa - Bak model without 
energy minimization (dataset SI), the shortest distances of Y89 Bak to 
R88 Bak , K35 Noxa , and R39 Noxa are 6.4, 3.7, and 7.6 A, respectively. These 
distances indicate cation-pi interactions of Y89 Bak with R88 Bak , 
K35 Noxa , and R39 Noxa , according to the survey of cation-pi interac- 
tions in protein crystal structures 30 . Therefore, the binding of Noxa to 
the Bak BH3-binding groove involves pi-pi interactions of F32 Noxa 
and cation-pi interactions of K35 Noxa and R39 Noxa with the respective 
Bak residues at the groove in addition to the anticipated hydrophobic 
and ionic interactions. 

In the BH3 domain of Bim, residues that correspond to F32 Noxa , 
K35 Noxa , and R39 Noxa are I155 Bim , E158 Bto \ and Y162 Bim respectively, 
raising the possibility that the pi-pi and cation-pi interactions are 
unique to the Noxa - Bak complex. Interestingly, the Bim - Bak model 
refined by microsecond-scale molecular dynamics simulations 
shows that F159 Bim , Y162 B,m , and Y163 Bim form pi-pi interactions 
with the Bak BH3-binding groove and extend the 10-residue aro- 
matic network of apo Bak to a 13-residue aromatic network in the 
Bim - Bak complex (datasets S4-S6). These residues also have cation- 
pi interactions with R88 Bak bridged by Y89 Bak (datasets S4-S6). 

In examining structures of anti-apoptotic proteins in complex 
with BH3-only proteins, we also found pi-pi and cation-pi inter- 
actions of residues at the BH3-binding groove of the following 
proteins: mPuma - mMcl-l with a 7-residue aromatic network 
(Protein Data Bank ID: 2ROC 21 ), mNoxa - mMcl-l with an 8-res- 
idue aromatic network (Protein Data Bank ID: 2ROD 21 ), mouse 
Bcl-X L in complex with mouse Bim possessing an 18-residue aro- 
matic network (Protein Data Bank ID: 1PQ1 31 ), and Norway rat 
Mcl-1 bound with Bim stabilized by a 9-residue aromatic network 
(Protein Data Bank ID: 2NL9 28 ). However, the pi-pi and cation-pi 
interactions are missing in Bfl-1 in complex with Noxa (Protein 
Data Bank ID: 3MQP) and mouse Mcl-1 bound to an alternative 
mouse Noxa BH3 domain (Protein Data Bank ID: 2JM6 28 ). These 
observations suggest that interactions of BH3-only proteins with 
the Bak BH3-binding groove as well as some of the anti-apoptotic 
Bcl-2 paralogs involve not only hydrophobic and ionic interac- 
tions but also pi-pi and cation-pi interactions. 
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Figure 3 | Main-chain conformational changes of Bak induced by Noxa 
binding, (a) Residues 51-61 blocking the groove between helices 1 and 6 in 
the unbound Bak crystal structure (Protein Data Bank code: 2IMT; Bak in 
gray), (b) Residues 51-61 partially vacating the groove between helices 1 
and 6 in the Noxa-Bak model obtained with energy minimization after 
molecular dynamics simulations (Bak and Noxa are in gray and green, 
respectively), (c) Overlay of the unbound Bak (magenta) and the Bak 
(gray) in complex with Noxa green) with arrows showing that Noxa 
slightly pushes helices a3 and a5 away, resulting in moving residues 51-61 
away from helix 5. 



Bak Homo-Oligomerization Induced by Bim or Noxa. Recently we 
reported the experimental observation that a high-affinity activator 
such as Noxa or Bim can bind Bak transiently to trigger Bak homo- 
oligomerization 14 . Knowing that the K D values of Noxa and Bim for 
binding to Bak in the presence of CHAPS are 24 nM and 29 nM, 
respectively 14 , one would expect hetero-oligomerization induced by 
these high-affinity binders. A homo-oligomerization mechanism at 
the structural level is, therefore, needed to resolve the paradox of tight 
but transient binding. 

Bcl-2 paralogs have been proposed to be part of a putative pore 
structure for cytochrome c release 32 34 . Three hypothetical homo- 
oligomerization mechanisms (asymmetric-single-conformer, asym- 
metric-two-conformer, and symmetric-single-conformer) have also 
been proposed for Bak 16,35,36 . To resolve the paradox described above, 
we investigated these mechanisms by building respective pore mod- 
els of Bak homo-oligomers. According to the symmetric-single-con- 
former mechanism, Bak everts its helix oc2 and subsequently forms a 
symmetric and reciprocal homo-dimer as a subunit of a Bak homo- 
oligomer 16 . However, we could not build a pore model for this mech- 
anism for two reasons. First, residues 51-71 of the first Bak clash with 
helices ocl, oc4, oc5, and oc6 of the second Bak, unless a large portion of 
helix ocl and/or helix oc2 of the first Bak is uncoiled. Second, the 
direction toward which the C-terminal transmembrane helix of 
one Bak points is the opposite of the other, which prevents half of 
the Bak transmembrane helices of the oligomer from being anchored 
in the membrane. We could not build a pore model for the asym- 
metric-two-conformer mechanism either. Alternatively, given the 
Bak conformational changes induced by Bim or Noxa, we built a 
Bak octamer pore model for the asymmetric-single-conformer 
mechanism that can resolve the paradox. 

To simplify the description of the Bak octamer model, we refer the 
two grooves between helices a3 and cx4 and between helices otl and 
a6 to as BH3 front -binding and BH3 back -binding grooves, respectively, 
because the groove between otl and oc6 is on the opposite side of Bak 
from the canonical BH3-binding groove. Formation of the Bak octa- 
mer presumably requires three steps (Figure 4a). Step 1, Noxa 
(indigo) binds an apo Bak (magenta) to partly expose the BH3 Back - 
binding groove of Bak 1 (red). Step 2, the resulting Noxa'Bak' com- 
plex binds another apo Bak to have the turn that joins helices oc2 and 
oc3 of Bak 11 (orange) to be docked at the BH3 back -binding groove of 
Bak 1 . This complexation involving the BH3 back -binding groove of 
Bak 1 partly exposes the BH3 back -binding groove of Bak 11 . Step 3, 
repeating Step 2 six times results in a hetero-nonamer (Noxa - 
Bak I -Bak II -Bak III -Bak IV -Bak v -Bak VI -Bak VII -Bak vm ) to partly expose 
the BH3 back -binding groove of Bak vm , which ejects Noxa and accom- 
modates the turn that joins helices oc2 and oc3 of Bak 1 to form a donut- 
shaped homo-octamer that is akin to the pore model mentioned in 
references 32-35 . 

Derived from the Noxa "Bak model with energy minimization and 
12 10-ns molecular dynamics simulations in explicit water, the Bak 
octamer pore model is supported by several experimental observa- 
tions. First, we have reported that the apparent molecular weight of 
the Bak homo-oligomer induced by Noxa is —200 kDa, ~8 times of 
the monomeric Bak molecular weight of 25 kDa 14 . Second, the Bak 
octamer is consistent with the report that detergent exposure induces 
formation of a Bax oligomer composed of 6-8 Bax monomers 33 . 
Although early indirect structural data 34,37-39 suggested that a Bax 
oligomer with ~4 or 7 Bax monomers might be responsible for the 
efflux of cytochrome C, a recent high-resolution (1.5- A) crystal 
structure of bovine cytochrome C 40 and studies using hydrodynamic 
analysis and atomic force microscopy 34 unequivocally show that the 
Stokes diameter of cytochrome C is much greater than the earlier 
reported diameter of 17 ± 3 A. Given the 90% sequence identity 
between human and bovine cytochrome C (National Center 
Biotechnology Information session numbers: NP_061820 and 
NP_001039526, respectively), the diameter of human cytochrome 
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Figure 4 | Proposed pore-forming Bak homo-oligomerization mechanism and its support, (a) A hypothetic pore-forming Bak homo-oligomerization 
scheme (the apo Bak is represented by a magenta sphere with the BH3-binding groove on the right and the groove between helices 1 and 6 on the left; the 
bound Bak is shown by a sphere in other colors; Noxa is depicted by a purple wedge), (b) An energy- minimized atomic model of Bak homo-octamer 
showing a cavity diameter suitable for the passage of the crystal structure of bovine cytochrome C (Protein Data Bank ID: 2B4Z; the heme ring is in stick), 
(c) Close-up view of three Bak units of the homo-octamer showing two interfaces of each unit of the Bak octamer that resemble the two interfaces of the 
Bax NMR structure (Protein Data Bank ID: 2K7W) with its C-terminal transmembrane domain (a9 in green) bound at the BH3-binding groove and a 
conformationally constrained peptide of the Bim BH3 domain (Bim SAHB in yellow) docked at a groove between two helices that correspond to helices al 
and a6 of Bak. 



C is at least 36.84 A according to the oval-shaped bovine cytochrome 
C crystal structure. A Bak heptamer has a calculated cavity diameter 
of —26 A that is too narrow for the passage of cytochrome C. The 
cavity diameter of the Bak octamer pore model is calculated to be 
~40 A (Figure 4b), hence indicating that the donut-shaped octamer 
has a cavity suitable for the passage of cytochrome C with a diameter 
in the range of 37-40 A and lending further credence to the Bak 
octamer model. 

Given the debate on whether the Bax/Bak-induced MOMP is 
selective 41 " 46 or nonselective 47,48 for effluxes of cytochrome C and 
other intermembrane-space proteins, and, more importantly, given 
the fact that other larger proteins than cytochrome C are also released 
from mitochondria during apoptosis, it is important to note the 
possibility that a series of —40- A holes punched by the Bak octamers 
aligned next to each other or aligned in a circle pattern can make a 
larger hole to release larger proteins such as Smac (second mitochon- 
dria-derived activator of caspases). 

Unlike the symmetric-single-conformer model, the octamer pore 
model permits all of the C-terminal transmembrane helices to point 



in the same direction and be orthogonal to the ring structure, thus 
enabling the ring structure to be attached to a membrane. Each Bak 
in the proposed octamer has two interfaces to make the pore struc- 
ture (Figure 4c), the BH3 fron '-binding groove at the interface on one 
side of each Bak and the BH3 back -binding groove at the interface on 
the opposite side. Consistent with this model, both front and back 
grooves (helices al, a3, 0(4, and tx6) are reportedly involved in Bak or 
Bax oligomerization 5,1415 . While our experimental studies showed 
that mutations of the BH3 fronl -binding groove can prevent Bak oli- 
gomerization 14 , others reported that the Bak activation also involves 
al 5 and oc6 49 and that the corresponding BH3 bilck -binding groove (al 
and a6) of a Bax NMR structure binds Bim SAHB , a stapled peptide that 
is thought to model the Bim BH3 domain, while the corresponding 
BH3 front -binding groove (a3 and a4) of Bax is occupied by its C- 
terminal transmembrane domain (a9) 15 . The two interfaces of the 
Bax NMR structure resemble the two interfaces of each Bak in the 
octamer (Figure 4c). This resemblance further strengthens the like- 
lihood that Bak homo-oligomerization involves a pore-forming Bak 
homo-octamer. 
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While computational refinement and further experimental veri- 
fication of this Bak octamer pore model are required, the proposed 
homo-oligomerization mechanism explains how high-affinity BH3- 
only activators can induce Bak homo-oligomerization. Because of the 
geometric restraint on both ends of the hetero-nonamer and the 
partial opening of the BH3 back -binding groove of Bak vm , the proximity 
and orientation effect 50 (viz., the entropic effect in a modern term) 
and the allosteric modulation greatly increase the affinity of the 
BH3 back -binding groove of Bak™ for the turn that joins helices oil 
and a3 of Bak 1 thus enabling the replacement of Noxa by the back 
groove of Bak™ 1 even though Noxa has a higher affinity for the 
BH3 front -binding groove of Bak 1 than apo Bak vm . In other words, 
when Bak 1 and Bak™ are restrained at both ends of the nonamer, 
the affinity of the BH3 back -binding groove of Bak™ for the turn of 
Bak 1 becomes higher than that of Noxa due to the partial opening of 
the back groove of Bak™ and the avoidance of the entropy loss for the 
complexation of Bak 1 with Bak™. This explains how a high-affinity 
activator such as Noxa or Bim can bind Bak transiently and trigger 
Bak foorao-oligomerization rather than hetero-oligomerization. 

In conclusion, simulation-refined 3D models of Noxa 'Bak and 
Bim* Bak show pi-pi and cation-pi interactions in addition to the 
anticipated hydrophobic and ionic interactions between BH3-only 
protein and Bak. The models also show that binding of Noxa or Bim 
to the BH3-binding groove exposes the groove between helices otl 
and oc6 that is on the opposite side of the Bak molecule from the BH3- 
binding groove. These observations, coupled with the reported 
experimental data, suggest that a BH3-only protein induces for- 
mation of a pore-forming Bak homo-octamer through tight but 
transient binding to the canonical BH3 -binding groove. 

Methods 

Model Preparation. All His, Glu, Asp, and Cys residues were treated as HIP, GLU, 
ASP, and CYS, respectively, for Bak in complex with Noxa or Bim and for the Bak 
octamer using the AMBER force field (ff99SB) 51 - 52 . The Noxa-Bak and BinrBak 
complexes were generated by manual docking as described above and by mutating 
Noxa to Bim, respectively. The Bak octamer was generated by manual alignment. The 
crystallographically determined water molecules were removed from the Bak crystal 
structure before the manual docking or alignment. The topology and coordinate files 
of the Noxa - Bak and Bim - Bak complexes were generated by the PREP, LINK, EDIT, 
and PARM modules of the AMBER 5 program {University of California, San 
Francisco). The topology and coordinate files of the Bak octamer were generated by 
the xLeap module of the AmberTools 1.5 program {University of California, San 
Francisco). The energy minimization was performed by using the SANDER modules 
of the AMBER programs Version 5 (for Noxa ■ Bak and Bim • Bak) and Version 1 1 (for 
the Bak octamer) with a dielectric constant of 1.0 and 500 cycles of steep est -descent 
minimization followed by 10,000 cycles of conjugate- gradient minimization. 

For Noxa - Bak, the docking- generated complex and the average complex structure 
of cluster 6 from the first-round simulations (Table SI) were used for the first- and 
second-round simulations, respectively; the energy- minimized complex was solvated 
by using EDIT with 5,897 and 6,744 TIP3P water molecules 53 for the first- and 
second-round simulations, leading to a system of 20,703 and 23,244 atoms, respect- 
ively. The mutation -derived Bim - Bak complex was solvated by using EDIT with 
6,599 TIP3P water molecules leading to a system of 22,773 atoms. The energy- 
minimized Bak octamer was solvated by using the xLeap module with 59,580 TIP3P 
water molecules leading to a system of 198,828 atoms. The water molecules were 
obtained from solvating the complex using a pre-equilibrated box of TIP3P mole- 
cules, whose hydrogen atom charge was set to 0.4170, where any water molecule was 
removed if it had an oxygen atom closer than 2.2 A to any solute atom or a hydrogen 
atom closer than 2.0 A to any solute atom, or if it was located further than 10.0 A along 
the x-, y-, or z-axis from any solute atom. 

Multiple Molecular Dynamics Simulations. The solvated system was energy- 
minimized for 100 cycles of steepest-descent minimization followed by 100 cycles of 
conjugate-gradient minimization to remove close van der Waals contacts in the 
system, then heated from 0 to 300 K at a rate of 10 K/ps under constant temperature 
and volume, and finally simulated independently with a unique seed number for 
initial velocities at 300 K under constant temperature and pressure using the 
PMEMD module of the AMBER program (University of California, San Francisco). 
For the AMBER program, Version 8 was used for Noxa* Bak and Bim - Bak, and 
Version 11 was for the Bak octamer. All simulations used (1) a dielectric constant of 
1.0, (2) the Berendsen coupling algorithm 54 , (3) a periodic boundary condition at a 
constant temperature of 300 K and a constant pressure of 1 atm with isotropic 
molecule-based scaling, (4) the Particle Mesh Ewald method to calculate long-range 
electrostatic interactions 55 , (5) a time step of 1.0 fs, (6) the SHAKE-bond-length 



constraints applied to all the bonds involving the H atom, (7) saving the image closest 
to the middle of the "primary box" to the restart and trajectory files, (8) formatted 
restart file, and (9) default values of all other inputs of the PMEMD module. All 
simulations were performed on a cluster of Apple Mac Pros with 92 Intel Xeon cores 
(3.0 GHz), a cluster of Apple Xserves with 590 G5 processors (2.2/2.4 GHz), and the 
Intel Xeon clusters at the University of Minnesota Supercomputing Institute. 

Data Analysis. Average structures were obtained by using the CARNAL module of 
AMBER 5. Cluster analyses were performed by using the PTRAJ module of 
AmberTools 1.5 (University of California, San Francisco). RMSDs and ACRMSDs 
were calculated by using the McLachlan algorithm 56 as implemented in ProFit V2.6 
(http://www.bioinf.org.uk/software/profit/) and the CARNAL module, respectively. 

For each of the 20 (round 1) or 52 (round 2) simulations of Noxa'Bak, 200 (round 
1) or 100 (round 2) instantaneous conformations were saved at 5-ps (round 1) or 10- 
ps (round 2) intervals during the last 1-ns period. A total of 4,000 (round 1) or 5,200 
(round 2) instantaneous conformations of Noxa - Bak from the 20 (round 1) or 52 
(round 2) simulations were subjected to a cluster analysis using the averagelinkage 
algorithm [epsilon = 3.0 A; RMS on Y89 Bak , F93 Bak , Flll Bak , F134 Bak , and F32 Noxa for 
Noxa-Bak or Y89 Bak , F93 Bak , Flll Bak , F134 Bak , and F159 Bim for BinrBak] 57 imple- 
mented in the PTRAJ module of the AmberTools package (Version 1.4 for Noxa 'Bak 
and Version 1.5 for BinrBak; University of California, San Francisco). Seven (round 
1) or 2 (round 2) clusters of the Noxa - Bak conformations were identified (Table SI). 
The average conformation of each cluster was energy- minimized using the SANDER 
module of the AMBER 5.0 program with a dielectric constant of 40.0 and 500 cycles of 
steepest-descent minimization followed by 10,000 cycles of conjugate-gradient 
minimization. The interaction energy between BH3-only protein and Bak of the 
energy- minimized complex was then calculated using an in-house progm with a 
dielectric constant of 40.0 and a nonbonded cutoff of 50,000 A (Table SI). 

Graphics. All figures were generated by using the MacPyMOL VI. 5.0 (Schrodinger 
LLC, Portland, OR), AutoCAD 2012 for Mac Version F.51.M.57 (Autodesk Inc, San 
Rafael, CA), and Adobe Photoshop CS5 Extended Version 12.1 x64 (Adobe Systems 
Incorporated, San Jose, CA) programs. 
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